RUB-TPII-6/96 
March 1996 
hep-ph/9604324 



Nucleon spin structure from the instanton vacuum^ 



C. Weiss^ 



Institut fiir Theoretische Physik 11 
Ruhr-Universitdt Bochum 
Bochum, Germany 



Abstract 

We discuss the evaluation of the nucleon isoscalar axial coupling, , in the in- 
stanton vacuum, using the 1/iVc expansion. This approach allows a fully consistent 
treatment of the U{\)a anomaly. We compute the nucleon matrix element of the 
topological charge density, {N\FF\N) ^ and show that it reduces to the matrix el- 
ement of the divergence of the isoscalar axial quark current. Our arguments show 
that the usual evaluation of in the chiral quark soliton model is consistent with 

the U{\)a anomaly in leading order of l/Nc- Such calculations give g^^ = 0.36, 
which is in agreement with the recent estimate by Ellis and Karliner. 
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The so-called "proton spin crisis" has attracted a lot of interest in recent years |1|]. 
The moments of the measured polarized proton structure function, gi, combined with 
information from neutron and hyperon beta decay, allows to extract the value of the 
isosinglet axial coupling constant of the nucleon, ^f^-*. This quantity is defined as the 
nucleon matrix element of the isosinglet axial current (we assume Nf quark flavors), 

J5,^{x) = 5I^/(3^)7/.75V'/(a;), (1) 
/ 

{N\J,,{0)\N) = gfui,i,u. (2) 

In the naive parton model, g^^^ can be interpreted as the fraction of the nucleon spin 
carried by the quarks, commonly denoted by AS. The value originally obtained by 
EMC was consistent with zero 0, while an analysis of newly available data by Ellis and 
Karliner comes to a value of g^^^ = 0.27 ± 0.04 ± . . . [Q. These results have prompted 
many theoretical investigations, which provided new insights about the relation of the 
parton model language to QCD. 

The question about a parton interpretation aside, it remains a challenge to understand 
the small value of g^^^ in QCD. Following the work of Brodsky et al. 0, who showed that 
9^A = in the Skyrme model, this quantity has been estimated in effective chiral models 
of the nucleon with explicit quark degees of freedom, such as the chiral bag model 0] and 
the chiral quark soliton model |^, which give non-zero results in subleading order of 
1/iVc. In these approaches one faces the question how the U{1) axial current operator 
of QCD, eq.(|l]), is to be represented consistently in terms of the degrees of freedom of 
the effective model. This problem does not arise for conserved currents related to "good" 
symmetries of the strong interactions, like the electromagnetic or isovector axial current. 
The f/(l) axial current of QCD, however, has an anomalous divergence, 

N ~ _ 
dMx) = Y^FFix) + 2i'£mf^n^^f{x), (3) 

which is non-zero in the chiral limit, mj ^ 0. Obviously, a consistent definition of the 
U{1) axial current in an effective description of QCD requires an understanding of the 
role of gluonic degrees of freedom — in particular, of topological fluctuations. 

In this note we want to show how a fully consistent and quantitative description of the 
nucleon isoscalar axial coupling can be achieved in the framework of the instanton vacuum. 
The relevance of instantons has been amply demonstrated by phenomenology as well as by 
lattice calculations 0, H. In particular, instantons describe the spontaneous breaking of 
chiral symmetry, which is the most important non-perturbative phenomenon determining 
the structure of light hadrons, including the nucleon p], |[. Using a l/A^^c-expansion, one 
derives from the instanton vacuum an effective quark action in the form of a Nambu- 
Jona-Lasinio model, which provides a convenient tool for computing correlation functions 
of meson and baryon currents. This leads to a description of baryons as chiral solitons 
— Nf. "valence" quarks moving in a self-consistent meson field |]10|. Such an approach 



gives a very good description of the static properties of non-strange and strange baryons 



as well as their various form factors [O 
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On the other hand, the instanton vacuum provides a microscopic picture of the non- 
perturbative fluctuations of the gluon field. It therefore offers the possibihty to include 
the effects of the U{1)a anomaly in an effective description of hadrons. In the meson 
sector, the solution of the so-called U{1)a_ problem — the large mass difference between 
the rj and rj' meson — has been one of the first successes of the instanton vacuum |T^ 



Recently, a method has been developed to evaluate nucleon matrix elements of gluon 
operators in the instanton vacuum. This method is based on the variational description 



of the instanton ensemble by Diakonov and Petrov |T3| and the 1/A^c^expansion in the 
quark sector. Within this approach, we can evaluate g^^^ in a way which is fully consistent 
with the U{1)a_ anomaly. The consistency may be demonstrated as follows. In QCD, the 
anomaly equation, eq.(^ is an operator equation, and g^^^ may equivalently be expressed 
as the nucleon matrix element of the topological charge density, FF[x), 

(iV|FF(0)|iV) = ^gfmr^uzj.u. (4) 



Using the method of we can now compute the nucleon matrix element of FF{x) 
in the instanton vacuum and show explicitly that the value of g^^ obtained from eq.(^) 
agrees with the nucleon matrix element of the usual isosinglet axial quark current. In 
other words, the instanton vacuum allows to evaluate the nucleon matrix element of both 
the left- and the right-hand side of the anomaly equation, eq.(|^), and we verify that both 
agree within our scheme of approximations. 



In the following, we give a brief outline of the approach developed in [14|. We first 
discuss how the [/(l)A^anomaly determines the dispersion of topological fiuctuations of 
the numbers of instantons and antiinsantons in the ensemble. We then formulate the 
prescription to calculate the nucleon matrix element of FF, and verify that the U{1)a 
anomaly is realized at the level of nucleon matrix elements. Finally, we consider the 
implications of this result for calculations of g^^ within the chiral soliton model [^, ^. 

U {1) A^anomaly and fluctuations of the numbers of instantons. In the instanton ap- 
proach the partition function of Euclidean Yang-Mills theory is reduced to a statistical 
ensemble of instantons (/) and antiinstantons (/). An essential property of this descrip- 
tion is that the number of "pseudoparticles" is not fixed (grand canonical ensemble). 
This is a necessary consequence of the trace and U{1)a anomalies of QCD. In fact, the 
anomalies unambiguously determine the dispersion of the fiuctuations of the number of 
/'s (A^+) and J's (A^-) in the ensemble. The dispersion of fiuctuations of the total nuber 



of instantons, A^^. + A^_, is governed by the trace anomaly (see |T^ for a discussion). 
Fluctuations of the difference of the number of /'s and J's, 

A = A^+-A^_ ^ 0, (5) 

are determined by the topological susceptibility of the vacuum, 

1 

32^2 

This follows from the fact that a single /(/) has topological charge ±1, 

^ ^ <fxFF{x)] _ = ±1, (7) 



(A^) = (Q?), Q, ^ I d'xFFix). (6) 



327r2 \J / single /(/) 
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if we assume that the VEV of is entirely due to instanton fluctuations. 

In QCD (with fermions), the topological susceptibility is in the chiral limit determined 
completely by the U{1)a anomaly. It is possible to derive the limiting behaviour by 
saturating the Ward identity for the correlation function of two topological charge densities 
with pseudoscalar meson states One finds {V is the four-dimensional volume) 

The topological susceptibility is proportional to the harmonic average of the quark masses 
and thus vanishes if one or more of the fermion flavors becomes massless. 

The instanton ensemble with fermions exhibits precisely this behavior. The vanishing 
of the topological susceptibility in the chiral limit is a consequence of the presence of 
"unbalanced" zero modes in the fermion determinant for 7^ A^_. In fact, a calculation 
of the fermion determinant in zero-mode approximation (keeping only the zero mode in 
the interaction of the fermions with instantons [Q) results in a probability distribution 
of A with a dispersion identical to eq.(^ Thus, the instanton vacuum correctly 

describes the topological susceptibility of the QCD vacuum in the chiral limit. 

Nucleon matrix element of FF. The instanton vacuum gives a quantitative prescrip- 
tion for evaluating the nucleon matrix elements not only of quark, but also of gluon 
operators. Specifically, we are interested in the zero-momentum transfer nucleon matrix 
element of the topological charge density, FF{x). It is extracted from the correlation 
function 

JNJiQt), (9) 



where = JAr(xi), j]^ = Jn{x2) are nucleon currents consisting of Nc quark fields 
coupled to a color singlet. The correlation function is to be computed as an average over 
the grand canonical ensemble of instantons with fermions. The averaging consists of two 
steps: 

I) average over the instanton coordinates and the fermion field of an ensemble with 
fixed N+,N^ (canonical ensemble), allowing for 7^ A^_. 

II) average over fluctuations of A = — A^_, with the dispersion (A^) = (Q'f) given 
by eq.(D. 

Step I: The fixed-A^± average can conveniently be carried out with the help of the 
effective quark action, which is obtained by integrating over the instanton coordinates 
first [1^, [1^ . It has the form of a Nambu- Jona-Lasinio modelQ, 

Nf 

S'eff[V^,V] = I d\Y.^f^xljf-Y+-Y^-S^, (10) 



^The fermion field used here is identical to —iil)^ of 
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where y± are 2A^/-fermion vertices of the form of the 't Hooft interaction (determinant 
in flavor), 



Nf-l 



y±[^M = (^) mff J d'xdetMx)f„ (11) 

J±i^)f9 = -^/^/^exp(-z(A;-/)-x)F(A;)F(/)V^;(/c)i(l± 75)^,(0. 

Eq. ([lO|) is derived by evaluating the fermion determinant of the ensemble of N± instantons 
in saddle-point approximation ( 1 /A^^c-expansion) . Here, M is the dynamical quark mass, 
and F{k) a form factor proportional to the wave function of the instanton zero mode. 
For unequal numbers of J's and /'s, the effective quark action contains a CP-violating 
violating fermion vertex proportional to A, 

This coupling of the quarks to topological fluctuations of the number of instantons, which 
emerges naturally in a saddle-point evaluation of the fermion determinant for A^_|_ 7^ N_, 
plays a crucial role in the realization of the U{1)a anomaly, as seen below. An effective 
action similar to eqs. ([lO| , |T2l) has also been suggested by Nowak et al. [|l^. These authors 
introduce a local (space-time dependent) density of instantons through a coarse-graining 
procedure. We do not require this concept in our approach. 

With the help of the effective quark action, eq.(|10|), the fixed-iV± average of two 
nucleon currents can be represented as 



VipVi) exp [-SeaW, j 

In the limit of large Euclidean times this correlation function is dominated by the saddle 
point corresponding to the static chiral soliton |TD|. Also correlation functions with gluon 
operators can systematically be expressed as averages in the effective quark theory |[T4 
In the case of the topological charge, eq.(^ the result is particularly simple, 

t n \ _ A / 7.. ft 



The free nucleon correlation function is simply multiplied by the topological charge of the 
fixed-A^± instanton ensemble. A, which seems intuitively plausible. 

Step II: Having computed the average over the fixed-A^-t ensemble, we now must 
average the result over fluctuations of the topological charge. A, 

JnJ^nQi) = , (15) 



. fixed-7V± 
A 



where -P(A) is the probability distribution of A in the grand ensemble, describing the 
dispersion eq.(^). Since the dispersion is 0{m) and thus small in the chiral limit, it is 
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sufficient to evaluate the fixed- average on the RHS of eq.(|T5|) to the first non-vanishing 
order in A. For our correlation function eg. ([llD this means 

= (A^)A(j,4)^^|^^^ (16) 
= (JnJKy^-Y^))^^ _ (17) 



A=0 



efT 



A=0 



Note that the factor (A^) has canceled with a corresponding factor in eg . ([T^) , so that the 
result is finite in the chiral limit. In eg.(|18|) we have used the fact the — is nothing 
but the divergence of the isosinglet axial current of the effective guark theory, eg. ([TO|) , 



d^xd,J,,{x) = 2Nf{Y+-Y_), (19) 

as follows from the eguations of motion in leading order of l/N^. 

Eg.(0) now shows how the ?7(l)A^anomaly is realized: the grand canonical correla- 
tion function with the topological charge reduces to the Nambu-Jona-Lasinio correlation 
function with an insertion of the divergence of the usual isosinglet axial current of the 
massive "constituent" guark field. Taking the limit of large Euclidean times, one obtains 
from eg.(|TBp the corresponding eguation for the nucleon matrix elements. 

Conclusions. We have demonstrated how the instanton vacuum naturally includes 
the U {1) Asmomalj in the effective description of the nucleon. The mechanism are the 
fluctuations of the difference of the number of J's and 7's in the ensemble, the dispersion 
of which is governed by the f/(l)A^anomaly. Our argument shows that a calculation of 
g^^^ in the chiral guark soliton model may be regarded as a consistent estimate of the 
nucleon matrix element of the topological charge, eg.(^), in the instanton vacuum. We 
have established this eguality in the leading order of the l/N^ expansion. A systematic 
treatment of higher-order corrections in l/N^. is in principle possible, but difficult. For 
instance, the 't Hooft form of the effective guark interaction, eg.(|l^), applies only to the 
leading order in l/Nc- 

Calculations in the chiral guark soliton model with SU{3) flavor give g^^^ ~ 0.36 ||^. 
This result is in agreement with the value obtained in the recent analysis by Ellis and 
Karliner 

The correct description of the U{1)a anomaly is a crucial test for the general method 
for evaluating matrix elements of gluon operators developed in |T^. This method allows 



also to evaluate the nucleon matrix elements of QCD operators of leading and non-leading 
twist, as are encountered in the OPE-description of deep-inelastic scattering. Such an 
approach opens the prospect of describing the deep-inelastic structure of the nucleon 
consistently with its hadronic structure, in one well-defined scheme of approximations. 
Calculations of the nucleon structure functions using this method are in progress. 



The ideas presented here have been developed together with D.I. Diakonov and M.V. 
Polyakov. It is my pleasure to thank them for a stimulating collaboration. I am also 
indebted to P.V. Pobylitsa for many interesting discussions. 
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